Isotropization in Brane Gas Cosmology 
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Brane Gas Cosmology (BGC) is an approach to unifying string theory and cosmology in 
which matter is described by a gas of strings and branes in a dilaton gravity background. The 
Universe is assumed to start out with all spatial dimensions compact and small. It has pre- 



' viously been shown that in this context, in the approximation of neglecting inhomogeneities 



and anisotropies, there is a dynamical mechanism which allows only three spatial dimensions 
to become large. However, previous studies do not lead to any conclusions concerning the 



o 

, isotropy or anisotropy of these three large spatial dimensions. Here, we generalize the equa- 

tions of BGC to the anisotropic case, and find that isotropization is a natural consequence 
of the dynamics. 



I. INTRODUCTION 



The Standard Big Bang (SBB) cosmology has become an extremely successful model that has been well 
tested by experiment. However, the model is incomplete. The underlying theory of classical general relativity 
^ i' and the description of matter as an ideal gas breaks down at the high temperatures of the early Universe, and 

the solutions of the theory in fact have an initial singularity. Moreover, SBB does not address many important 
cosmological questions such as the observed homogeneity, spatial flatness, and the origin of structure in the 
Universe. Cosmological inflation (see e.g. [1,2] for textbook treatises and [3,4] for shorter reviews) builds 
on SBB cosmology providing a solution to some of these issues, but it (at least in the context of scalar 
field-driven inflation) suffers from the same initial singularity problem [5] and other conceptual problems [4] , 
which indicate that inflation cannot be the complete story of early Universe cosmology. 

In recent years, many models motivated by string theory and M-theory have emerged as possible solutions 
to the outstanding problems of early Universe and inflationary cosmology (see e.g. [6,8] for recent but 
incomplete reviews). Beginning with the work on Pre-Big-Bang cosmology [9,10] it was realized that a 
dynamical dilaton should play an important role in the very early Universe. More recently, models have 
become prominent in which our Universe consists of a 3-brane embedded in a higher dimensional bulk space, 
with the standard model constrained to live on the brane [11-17]. Although these models can resolve a 
number of issues, such as the hierarchy problem, they introduce several other difficulties in the process. 
For example, large extra dimensions should be explained by classical general relativity, and it has been 
shown this results in problems stabilizing the brane [18]. More importantly, in most of these models the six 
"extra" spatial dimensions are taken to be compactified, a priori, with no explanation provided for how this 
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could come about dynamically. Although it appears to be an important concern for the naturalness of these 
models, this issue is rarely discussed in the literature. 

An alternative approach to string/M-theory cosmology is the string gas or BV scenario. This model 
began with works [19,20] in which the effects of string gases on the cosmological evolution of the low energy 
effective string theory background geometry including the dilaton were explored. The most important result 
to emerge from these works is a dynamical mechanism, tied to the existence of string winding modes, which 
yields a nonsingular cosmology and may explain why at most three spatial dimensions can become large if 
the initial state is chosen to correspond to a Universe which is small in all spatial directions. 

This work has been generalized to include the cosmological effects of p-brane gases and leads to the current 
model of Brane Gas Cosmology (BGC) [21]. In BGC, the Universe starts analogous to the SBB picture, i.e. 
hot, dense, and with all fundamental degrees of freedom in thermal equilibrium. The Universe is assumed to 
be toroidal in all nine spatial dimensions and filled with a p-brane gas. The assumption of toroidal geometry 
of the background space leads to the existence of string winding modes, since the background space admits 
cycles on which branes of the relevant dimensionalities (in particular one branes) can wrap. This wrapping is 
associated with a winding energy which - in the context of dilaton gravity - acts as a confining potential for 
the scale factor preventing further expansion of the spatial dimensions. Also associated with the brane arc 
oscillatory modes described by scalar fields and momentum modes which correspond to the center of mass 
motion of the brane. The momentum modes are related by T-duality to the winding modes and this duality 
results in the non-singular behavior of the model. In order for dimensions to decompactify, p-brane winding 
modes must annihilate with anti-winding modes and it is argued that this only occurs in a maximum of 
2p+l dimensions [21]. Since strings {p = 1) are the lowest dimensional objects that admit winding modes, 
since they are the lightest of all winding modes and hence fall out of equilibrium later than other winding 
modes, it follows that the number of large space-time dimensions can be at most (3+1). 

In the context of the background equations of dilaton gravity, the winding modes yield a confining potential 
for the scale factor which also gives rise to a period of cosmological loitering (expansion rate near zero) for 
the three large spatial dimensions. This is due to the time needed for winding modes to annihilate and 
produce closed strings or loops [22]. This is of great interest since loitering can explain the horizon and 
relic problems of standard cosmology without resorting to inflation ^. It was also shown in [22] that by 
considering loop production at late times BGC naturally evolves into the SBB, with a 3 -|- 1 dimensional, 
radiation dominated Universe. 

There are important issues that remain to be addressed within BGC. The fact that toroidal geometry was 
assumed for the background space was used for the existence and topological stability of winding modes. 
However, K3 or Calabi-Yau manifolds are more realistic choices for backgrounds within string theory and 
they do not admit 1-cyclcs (necessary to have topologically stable winding modes). Promising results have 
recently appeared which indicate that the conclusions of BGC extend to a much wider class of spatial 
background, including backgrounds which are K3 or Calabi-Yau manifolds [7,24]. 

Perhaps the main issue to be addressed in BGC is that of spatial inhomogeneities. That is, we would 



'^However, to obtain a solution of the flatness and entropy problems, a phase of inflation following the decoupling 
of the three large spatial dimensions may be required 
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naturally expect fluxes and p-brane sources to lead to the possibility of catastrophic instabilities of spatial 
fluctuation modes. Although we do not address this issue here, we plan to study the role of inhomogeneities 
in foUowup work. 

Other important issues for BGC arc the questions of stabilization of the six small extra dimensions and 
isotropization of the three dimensions that grow large. Although these topics may seem to be unrelated, they 
both can be addressed by generalizing the BGC scenario to the anisotropic case. This paper will concentrate 
on the isotropization of the three large dimensions, but the generalization of BGC to the anisotropic case 
achieved in this paper will be valuable to address the issue of stabilization in later work. 



We begin with a brief review of BGC, for more details the reader is referred to [21]. Consider compact- 
ification of 11-dimensional M- theory on S^, which yields type II- A string theory in 10-dimensions. The 
fundamental degrees of freedom in the theory are 0-branes, strings, 2-branes, 5-branes, 6-branes, and 8- 
branes. The low energy effective action of the theory is that of supersymmetrized dilaton gravity. 



where G is the determinant of the background metric 0^,^, (j) is the dilaton given by the radius of the 
compactification, H denotes the field strength corresponding to the bulk antisymmetric tensor field B^,^, 
and K is determined by the ten-dimensional Newton constant. 

Fluctuations of each of the p-branes are described by the Dirac-Born-Infeld (DBI) action [25] and are 
coupled to the ten-dimensional action via delta function sources. The DBI action is 



where Tp is the tension of the brane, gmn is the induced metric on the brane, bmn is the induced antisymmetric 
tensor field, and Fmn is the field strength tensor of gauge fields living on the brane. The constant a' ~ 1"^^ 
is given by the string length scale 1st- 

In this paper we will concentrate on the role of fundamental strings {p = 1), since strings play the critical 
role in the dynamics [21]. We will also ignore the effects of fluxes and string oscillatory modes since we 
are working in the low energy regime. Our basic approach is to consider the string gas as a matter soiirce 
for the dilaton-gravity background. The string winding modes play a key role. The main point is that we 
want to find the first departures from the standard FRW picture resulting from considering stringy effects. 
This should result from the winding modes of the strings and the presence of the dilaton in the equations of 
motion. 

In this paper, we generalize the equations of BGC to the anisotropic case. Thus, we take the metric to be 
of the form 



II. BRANE GAS COSMOLOGY 






(2.2) 



D 




(2.3) 



i=l 
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where the label i runs over the D spatial indices, Xi are the comoving spatial coordinates, t is physical time, 
and the scale factor in the i'th direction is log{Xi). We again stress that inhomogenities are of vital impor- 
tance, however we will leave their investigation to future work and take the string gas to be homogeneous. 
By varying the total action we obtain the following equations, 

D 



-^A, + ^2^e'^£ (2.4) 
-Y^X, + ^= e^E (2.5) 



=1 



Xi - ifXi = ^e^Pi, (2.6) 

where E is the total energy of the strings and Pj is the pressure in the ith direction, and we have introduced 
a "shifted" dilaton field 

D 

<^ = 2</,-^Ai. (2.7) 

i=l 

We can see immediately from the first equation that (p can not change sign. To insure that the low energy 
approximation remains valid wc choose ip < and (p > 0. This causes the dilaton to have a damping effect 
in the equations of motion ^ . Another important observation is that these equations reduce to the standard 
FRW equations for fixed dilaton as was discussed in [20]. This is comforting because at late times it is 
expected that the dilaton will gain a mass, perhaps associated with supersymmetry breaking and as a result 
we naturally regain the usual SBB cosmology. The string winding states will appear at late times as solitons 
[21]. 



A. Winding Modes 



The energy and pressure terms follow from varying the brane action (2.2). The energy associated with 
the winding modes in the ith direction is given by (taking the winding strings to be straight for simplicity) 

= nNi{t)ai{t) = fiNi{t)e^% (2.8) 

where A^, is the number of winding modes in the ith direction and /x is the mass per unit length of the string 
multiplied by the initial spatial dimension. It follows that the total energy is given by, 

D 

E}^=Y,P^Ni{t)e^K (2.9) 

i=l 



^Alternatively we could consider > and one might anticipate that this could lead to inflation. However this is a 
false conclusion, since the presence of winding modes act as a confining potential which can be seen from the negative 
pressure term in equation (2.6). 



4 



In the same approximation that the winding strings are straight, the corresponding pressure terms are given 

by [21], 

P^^ = -fiN,{t)e^^ (2.10) 
Pr = • (2.11) 

It follows then by inserting (2.10) into (2.6) that winding strings give rise to a confining potential for the 
scale factor, and hence prevent expansion. In order for a spatial dimension to be able to become large, the 
winding modes in that direction must be able to annihilate. 



B. Annihilation and Loop Production 

As argued earlier, the fact that strings are likely to annihilate in a maximum of three space dimensions 
leads to three of the initial nine spatial dimensions growing large. As the background continues to expand 

the three dimensional space will be filled with loops resulting from the annihilation of the winding states. 
The corresponding energy loss in winding states and the energy and pressure of the creation of loops must 
be taken into account when considering the dynamics. The strings in the expanding space will become of 
macroscopic size and the required pressure and energy terms will be analogous to that of a cosmic string 
network [26,27]. 

We know that the energy in winding modes is given by (2.9). Considering the time rate of change of this 
energy we find, 

D D 

El^ = ^lxNi{t)e^' +J2iiNi{t)e^^Xi. (2.12) 

1=1 i=l 

The first term on the right hand side of this equation corresponds to energy loss into string loops or radiation, 
the second term corresponds to the gain in total energy of the strings due to the stretching by the expansion 

of space. 

The energy in loops can be written as [22] 

E^oops ^ g^^^y^^ (2.13) 

where g{t) is the energy density in string loops (plus radiation) per initial comoving volume (and is a constant 
if no loop production or energy loss from winding strings into radiation occurs), and Vq = exp(^Aoi) is 
the initial volume, Aoi being the initial values of the logarithms of the scale factors. Since sting loops are 
likely to be produced relativistically, we will use the equation of state of radiation to describe them. This 
approximation also allows us to treat the string loops and other radiation together. This approximation could 
easily be relaxed without changing our basic conclusions. Thus, we use the equation of state of relativistic 
radiation p = {l/3)p (p and p denoting energy density and pressure, respectively) which implies, 

ploops ^ ^g^t)V0. (2.14) 

To find how g{t) evolves we equate the loss in winding energy due to energy transfer (the first term on 
the right hand side of 2.12) with the change in loop energy, 
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D 

-Y,I^Ni{t)e^' =g{t)Vo- (2-15) 

i=l 

Thus, the loop production is determined by the change in the number of winding modes, as expected: 

D 

3(f) = -y-i^^iV.(f)e^*. (2.16) 

i=l 

It remains to determine the rate of winding mode annihilation, but before doing so let us consider the specific 
case of a 2+1 anisotropic Universe. 

III. ANISOTROPIC GENERALIZATION 

As three of the nine spatial dimensions grow large there is no a priori reason to expect that this should 
happen in an isotropic manner. Moreover, since the dimensions expanding correspond to the annihilation of 
winding modes due to string intersections, we might expect this process to occur at differing rates leading 
to anisotropic dimensions. However, we might also expect that as the winding modes annihilate in one 
dimension and that dimension begins to expand faster there arc then fewer winding modes left to annihilate. 
In this way the remaining dimensions are given an opportunity to isotropize. To explore if this is indeed the 
case, let us consider the case of a 2+1 anisotropic background, whose metric we write in the form 

ds^ = - e^\dx^ + dy^) - e^^dz^ 

The scale factor corresponding to A is denoted a{t), the one corresponding to v is denoted h{t). 
The equations of motion (2.4) then become, 

-2>? - ?>2 + = e*'^ 
-2A2-!>2 + ^ = ie^^ 

A - = ^-e^Px 

The energy and pressure terms are given by 

E = E:^ + E^°°p' = 2nN{t)e^ + fiM{t)e'' + g{t)e^^°+''° , (3.6) 
Pa = -l^Nit)e^ + ^g{t)e'^o+u° , (3.7) 

P. = -fiM{t)e'' + ^git)e^^"+-'' , (3.8) 

where N(t) and M(t) are the number of winding modes in the A and u directions, respectively. 

We now consider the effect of loop production as a result of winding mode annihilation. To simplify the 
analysis, we will assume that winding modes in A and v directions form two separate, noninteracting gases. 
We expect that this approximation will reduce the isotropization of winding modes, and hence that including 
the interactions we omit will lead to accelerated isotropization. This issue is being studied currently [28]. 
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(3.1) 



(3.2) 
(3.3) 

(3.4) 

(3.5) 



Winding mode annihilation results from string intersections, which depends [26,27] on the square of the 
number of modes, inversely on the cross-sectional area available for the interaction and is proportional to 
the Hubble length of the string. Since causality plays a role in the interaction, we must scale our dimensions 
by the Hubble time t. Hence, 

^(^^ = -'^""H^y = -'-''H^y = --^^^We-^-^ . (3-9) 

where cjv is a dimensionless constant, and cjv is this same constant rescaled by the basic length dimensions 
related to a and b. Similarly, for M{t) we have 

M(t) = -CMtM^(t)e-^^. (3.10) 
Using these expressions in (2.16) gives us the evolution of loop production, 

<?(t) =/xte-'^°-"(2cjviV2(t)e-" + CMM2(t)e-2^+"). (3.11) 



For computational simplicity we define l{t) = A, q{t) 
set of first order ordinary differential equations. 



= !>, and f{t) = (p, which leaves us with the following 



(3. 


,12) 


(3. 


,13) 


(3, 


,14) 


(3. 


,15) 


(3, 


,16) 


(3, 


,17) 


(3. 


,18) 


(3. 


,19) 


(3. 


,20) 


(3. 


,21) 



f = 2f+q^ + ie'^ (2fiN{t)e^ + nM{t)e'' + g{t)e^^<'+''°) 

q = W+ \e^[-I^M(t)e^ + ^-g{t)e^>^»+-^>) 
Nit) = -CNtN'^{t)e-^-'' 
M{t) = -CMtM^{t)e-'^^ 
g{t) = tite-^^"-" (2cNN^{t)e~'' + CMM^{t)e~^^+'')v^^ 

l{t) = A 
q{t) = 
fit) = 'p. 



The system (3.12-3.21) is over-determined (there are more equations than unknowns). Therefore we take 
(3.12) as a constraint on the initial conditions. Wc find from numerical analysis that the winding modes 
N(t) and M(t) approach zero as the system evolves in time, as expected from winding mode annihilation 
(see Fig. 1). As a result loop production continues until all the winding modes are annihilated at which 
time g{t) approaches a constant (see Fig. 2). Note from Figure 1 that before converging to zero, the 
two winding numbers converge to eachother. This is a necessary condition for isotropization. This result 
seems reasonable, since the rate of annihilation depends on the number of winding states present. Thus, as 
the number of winding modes decreases in one direction, the other direction should have an larger rate of 
annihilation. This suggests a sort of equilibration process due to the presence of winding modes. 
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FIG. 1: A plot of the numbers A'^ (upper 
curve for small values of t) and M (lower 
curve) as a function of time t. We see that the 
number of winding modes approaches zero as 
winding annihilation continues. 




FIG. 2: A plot of the comoving energy den- 
sity in loops g as a function of time t. As 
the winding modes annihilate, loop produc- 
tion continues until all the winding modes 
have vanished and the loop energy becomes 
constant. 
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FIG. 3: A plot of the two expansion rates 
I (upper curve for small values of t) and q 
as a function of time t. It is seen that as 
the winding modes annihilate, the expansion 
rates converge to each other before converg- 
ing to zero. 
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FIG. 4: A plot of the anisotropy parame- 
ter ^ as a function of time t for a range of 
initial anisotropics. Curve a represents no 
initial anisotropy, but varying winding num- 
bers. Curve h has an initial anisotropy of 
A = 0.025 and curve c has y4 = 0.160. We 
see that the anisotropy parameter reaches a 
maximum early on in the evolution and in all 
cases tends to zero at late times. 




However, the above alone is not enough evidence to prove that isotropization of the dimensions occurs. To 
prove that isotropization occurs we need to study the geometry of the background. We can see from Figure 
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3 that the expansion rates, l{t) and q{t), converge at late times. Note, however, that at the same time the 
ratio of scale factors continues to increase. A more quantitative definition of isotropization can be obtained 
by defining the average Hubble parameter A and the anisotropy parameter A [29] , 

1 ^ 

A^-^Ai (3.22) 
-^-^E^^- (3-23) 

i=l 

In our case we have D = 3 and the anisotropy parameter becomes, 

^-\m+q{t)r ^ 

We find that for any amount of initial anisotropy, A approaches a maximum value and then goes to zero, as 
can be seen from Figure 4. Furthermore, for the case of isotropic initial expansion but inequivalent winding 
mode numbers we find that A reaches a larger maximum but the conclusion of isotropization at later times 
is unchanged. 



IV. CONCLUSION 

We have generalized the equations of BGC to the anisotropic case including the effects of winding state 
annihilation and loop production. We address the issue of isotropization of the three large dimensions 
quantitatively by introducing the anisotropy parameter A. Our analysis indicates that for an arbitrary 

amount of initial anisotropy, the anisotropy will reach a maximum early in the evolution and then approach 
zero at later times. Thus, we have explained how isotropy arises as a natural consequence of BGC. 
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